Efficient ion acceleration by collective laser-driven electron dynamics with ultra-thin 

foil targets 
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Experiments on ion acceleration by irradiation of ultra-thin diamond-like carbon (DLC) foils, with 
thicknesses well below the skin depth, irradiated with laser pulses of ultra-high contrast and linear 
polarization, are presented. A maximum energy of 13 MeV for protons and 71 MeV for carbon ions is 
observed with a conversion efficiency of ~ 10%. Two-dimensional particle-in-cell (PIC) simulations 
reveal that the increase in ion energies can be attributed to a dominantly collective rather than 
thermal motion of the foil electrons, when the target becomes transparent for the incident laser 
pulse. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.50.Jm, 52.65.Rr 
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Recent experiments in the field of relativistic laser- 
plasma interaction have shown that the conversion ef- 
ficiency (CE) from the laser to the kinetic energy of the 
ion bunch as well as the maximal energy of the ions can 
be improved by the use of ultra-thin foil targets in the 
range of 100 nm 

[EH llll- Additionally, it was empha- 
sized that besides the standard target normal sheath ac- 
celeration (TNSA) other, radiation pressure domi- 
nated acceleration mechanisms become possible for ultra- 
thin targets 0, B 0- An ultra-high contrast is required 
to avoid substantial expansion of the targets before the 
interaction with the main pulse. These conditions can 
be achieved by the use of a specially designed laser sys- 
tem in combination with double plasma mirrors (DPM) 
[E E) II | • The partial transmission of the intense laser 
pulse through an expanding target is expected to play 
a decisive role QJi, [HI [lj] and has been lately demon- 
strated in the regime of relativistic transparency [l5| . 

In this letter we present experimental results on ion 
acceleration from DLC foils of thicknesses ranging from 
50 nm down to 2.9 nm. The targets are irradiated by lin- 
ear polarized pulses of 45 fs FWHM duration focussed 
to a peak intensity of up to 5 x 10 19 W/cm 2 . We find 
an optimum in ion acceleration at a target thickness of 
5.6 nm, where ion energies reach values of 13 MeV for 
protons and 71 MeV for carbon ions. Thus, for the ultra- 
short pulses used, the optimum foil thickness is well be- 
low the collisional skin depth l s ~ 13 nm of the heated 
target, which is thus becoming transparent to the laser. 
The corresponding CE is reaching values of ~ 1.6% for 
protons (E > 2MeV) and - 10% in the case of C 6+ 
(E > 5MeV). By means of 2D PIC-simulations we find 
that the transmitted laser field imposes a dominantly col- 
lective rather than thermal motion of the foil electrons, 



leading to increased ion energies and enhanced CE. Our 
findings are supported by a semi-analytical model, which 
shows good agreement with the experimental results. 

The experiments were performed at the MBI - TW 
Thsapph laser of central wavelength 810 nm delivering 
1.2 J in 45 fs FWHM pulses with an amplified sponta- 
neous emission (ASE) contrast ratio smaller than 10 -7 up 
to ~ 10 ps prior to the arrival of the main peak. By means 
of a re-collimating DPM [l| , this contrast was increased 



by estimated four orders of magnitude [ll| , which is es- 
sential for the suppression of pre-heating and expansion 
due to the pulse background. The energy throughput of 
this DPM system was improved to values of 60-65 % com- 
pared to our previous experiments [H, resulting in pulse 
energies of 0.7 J. Finally, the laser pulse was focused on 
the DLC target with a f/2.5 parabolic mirror down to 
6 fim diameter and is diffraction limited by 30% under 
normal incidence. This corresponds to a peak intensity 
of 2.6 x 10 19 W/cm 2 or a normalized vector potential of 
ar, = 3.6. For a second set of experiments the focussing 
procedure was improved, leading to a focus diameter of 
3.6 fim and therefore a peak intensity of 5 x 10 19 W/cm 2 
or at = 5. The resulting ions were reg istered with a 
MCP coupled to a Thomson-Parabola [lq |. 

DLC targets of thicknesses ranging from 2.9 — 50 nm 
were used, having a density of 2.7 g/cm 3 . Owing to the 
high fraction of sp 3 -, i.e. diamond-like bonds of ~75%, 
DLC offers unique properties for the production of me- 
chanically stable, ultra-thin, free standing targets, such 
as exceptionally high tensile strength, hardness and heat 
resistance. The thickness of the DLC foils was character- 
ized by means of an atomic force microscope (AFM), in- 
cluding the hydrocarbon contamination layer on the tar- 
get surface which was present during the experiments. In 
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FIG. 1: Maximum energy per atomic mass unit for both pro- 
tons and carbon ions plotted vs. foil thickness for values of the 
normalized laser vector potential of arj = 3.6 and arj = 5. The 
experimental results are in excellent agreement with numbers 
deduced from 2D-PIC-simulations. 



addition, in order to precisely determine the structure of 
the contamination layer, the depth-dependend composi- 
tion of the foil was measured via Elastic Recoil Detection 
Analysis (ERDA). From these measurements we obtain 
a thickness of ~ 1 nm for the hydrocarbon contamination 
layer. Throughout the manuscript we are referring to the 
combined thickness of bulk and surface layer as it appears 
in the actual ion acceleration experiment presented. 

The maximum detectable energy/nucleon (E max ) ob- 
tained for protons and fully ionized carbon ions C 6+ is 
plotted as a function of the target thickness in FIG. Q] 
for two different laser intensities. The cut-off energies for 
both ion species exhibit a strong dependence on target 
thickness. In case of protons and ao = 5, E m:ix increases 
from around 7.5 MeV for a 40 nm foil up to 13 MeV for a 
5.6 nm foil, while for C 6+ the maximum energy rises from 
26 MeV for use of a 40 nm foils to 71 MeV for 5.6 nm. 
Further decrease of the target thickness down to 2.9 nm 
results in a steep drop of the observed ion energies. The 
energy distributions of all species are continuous as ex- 
pected. However, especially the observed carbon C 6+ 
energy of 71 MeV reaches for the first time a range of 
values that were previously only accessible by large sin- 
gle shot Nd:glass laser systems with 30 — 50 J pulse energy 
lil 17 1 . The shot-to-shot variations of 10%, indicated 
by the error bars, arise mainly from fluctuations of the 
laser pulse itself and from macroscopic modulations of 
the target surface. 

Furthermore, the CE (FIG[2]) was calculated by nu- 
merically convoluting an energy dependent divergence 
of the accelerated ion beam with the initially measured 
(ao = 5) spectra. These values (cf. inset in FIG(2]) were 
extracted from PIC simulations and supported by exper- 
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FIG. 2: Calculated CE for protons and C 6+ ions as a function 
of the target thickness, based on a energy dependent diver- 
gence angle of the ion beam, which was extracted from PIC 
simulations and supported by measured beam profiles (cf. fig- 
ure inset). The experimental data are in good agreement with 
the values of the CE obtained from PIC simulations. 



imentally obtained proton beam profiles using a stack of 
radiochromic film (RCF) layers. Comparing those to the 
divergencies of ion beams generated by other laser sys- 
tems e.g. 



18(, our divergence is smaller and the linear 
dependence on the energy is reasonable in the considered 
high-energy part. The dependence of the CE on the tar- 
get thickness is showing analogous characteristics to E max 
(cf. FIG.l). A sudden rise is observed when entering the 
laser transparency regime (D < l s ) reaching 10.5% for 
C 6+ and 1.6 % for protons at the optimum thickness. The 
CE drops down to below 1% when the target thickness 
reaches the TNSA dominated regime 0, H3| ■ 

The experiments were compared to 2D-PIC simula- 
tions, where the laser pulse was modeled by a Gaussian 
shape in time with a FWHM of 16 laser cycles, a Gaus- 
sian intensity distribution in focus with a FWHM spot 
size of 4 /xm and an ao = 5 (ao = 3.6). It interacts with a 
rectangular shaped plasma of initial density n e = 500 x n c 
consisting of 90 % carbon ions and 10 % protons (in num- 
ber density) to account for the presence of a contamina- 
tion layer. The simulation box is composed of 1200 x 
10000 cells with 1000 particles per cell and a total size 
of (10 x 20) /urn 2 . The total simulation time r given in 
laser cycles is 120. 

FIG. 1 shows that the simulated maximum proton 
energies as well as the carbon energies are in excel- 
lent agreement with the experimental data. In partic- 
ular, the optimum target thickness of 5.6 nm is repro- 
duced. This thickness for the peak ion energy is consis- 
tent with the thickness given by the empirical relation 
! 

c ~ 3 + 0.4 x op = 5 that was found in multiparametric 
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PIC-simulation studies by Esirkepov et al. [lj|. Here, a 
normalized areal electron density a — (n e /n c ) ■ (D/\l) is 
introduced, where D is the target thickness, A^ the laser 
wave length, n e the electron density of the target and 
n c the critical electron density, which in our case calcu- 
lates to a — 4.6. For foil thicknesses below the optimum 
(<t < 3 + 0.4 x a ), the plasma becomes increasingly trans- 
parent and the pulse is more transmitted than absorbed. 
Due to the low number of electrons in the focal volume 
(~ 10 11 ) their electric current is no longer sufficient to 
(i) reflect the laser pulse and (ii) to establish an effec- 
tive longitudinal charge separation field. This results in 
a sudden drop in ion energies and 50 % of the CE (the 
same amount as the reduction of the target thickness), 
as it was observed in the experiment (see FIG. 1+ 2). In 
case of a > 3 + 0.4 x a , the laser intensity is not suf- 
ficient to generate the maximum possible displacement 
of all electrons within the focal volume which gives rise 
to a decrease in the longitudinal charge separation field. 
Note that the optimum thickness for ion acceleration pre- 
dicted theoretically and observed experimentally is much 
smaller than previously used target thicknesses. 

We compared our results to a number of analytical 
scaling laws that have been derived for the standard sce- 
nario of ion acceleration from the back surface of a thin, 
but opa que foil, i.e. target normal sheath acceleration 
(TNSA) [I], EH, El. However, we find that all of the ex- 
isting analytical models fail to even qualitatively predict 
the two main features seen experimentally, namely the 
distinct peak in ion energies at a » 3 + 0.4 x ao as well 
as the strong dependence of the cut-off value on target 
thickness. Consequently, a refined theory is needed. 

In FIG. 3a the observed carbon energy spectrum is 
shown along with the simulated one agreeing with each 
other in their overall spectral shape. To illustrate the un- 
derlying ion acceleration mechanism in our case we show 
in FIG.[3];,d detailed electron and carbon densities after 
40 laser cycles as a function of the laser propagation di- 
rection x and the transverse coordinate y in units of the 
laser wave length Al . It is immediately striking that ions 
are accelerated strongly asymmetric heavily favoring the 
direction of laser propagation. This is in strong contrast 
to the model based on the self-consistent solution of the 
Poisson equation presented in [22j which predicts a sym- 
metric acceleration that happens primarily after the end 
of the laser pulse. From FIG. 3c+d it can be seen, that 
the carbon ions are accompanied by co-moving electrons 
which accelerate the ions in forward direction. 

Unlike the TNSA consideration in case of thicker tar- 
gets, in our current case, coherence of laser-driven elec- 
tron motion manifests importantly. In the typical TNSA 
with thicker targets, none of the laser field is transmit- 
ted as the foil remains opaque throughout the interaction 
duration. The electron motion is an indirect result of in- 
teraction with the laser at the front surface of the target 
and subsequent multiple interactions with matter in the 
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FIG. 3: (a) Comparison between experimentally obtained and 
simulated C 6+ spectrum as deduced from a 2D-PIC simula- 
tion of a 4.5 nm foil target (b) Electron propagation angle 
(blue circles) and laser transverse electric field E y (red line) 
at time t= 30 r, plotted vs. laser propagation direction x, 
given in units of the laser wavelength A. r is given in laser 
cycles. The green line at x — 5A gives the initial location of 
the target. For the same simulation, (c) gives the electron 
and (d) the carbon density distribution vs. x at time t= 40t. 



target bulk. That scenario leads to thermal motion as the 
electrons traverse the target to set up a quasi-static field 
at the back. In contrast, in nm-scale foils of thickness be- 
low the skin depth the laser imprints a directed motion 
onto electrons as those targets are partially transparent. 
These coherent electron motions can not be treated as 
thermal motion, but rather as collective motion tied to 
the laser and thus directly to the laser amplitude ao- The 
propagation angle of electrons upon laser interaction with 
the target as extracted from the 2D PIC simulation at 
time 30 r is shown in FIG. 3b. Obviously, the electron 
dynamics are classified in three categories, (i) electrons 
are driven forward by laser, (ii) electrons are forced to 
return backwards due to the electrostatic field, (iii) elec- 
trons are performing a collective, transversal oscillation 
owing to the laser-created coherent wave structure. Ad- 
ditionally, the laser transmission coefficient can be de- 
termined to be T ~ 0.1, which agrees with a theoretical 
value of T = 1/[1 + (ir<r/j p ) 2 ], where j p is the Lorentz 
factor in the laser \fa? L + 1 [23| . 

Initially, without assuming thermalized electrons (i.e. 
in contrast to j24[), we now formulate the maximal ion 
energies obtained in the laser driven foil interaction in a 
semi-analytical model based on [25j |. 

Following the analysis of Mako & Tajima [23] and 
Yan et al. (26[, the forward current density of elec- 
trons J and electron density n e are related through 
J = -e J^ mm V x gdV x and n e = § J max ^{^dv, where 
g is the electron distribution function and V max = 
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Here we note that at a given 



position the electron energy Eq in the electrostatic po- 
tential (j> is given as an implicit function of space and 
time through e = 7m e c 2 — ecj). We find that J(e) = 
—Jo (1 — e/Eo) a where a, a measure for the coherence is 
found to be about 3 in our simulation for a target thick- 
ness of 4.5 nm. Considering the relativistic dynamics of 
electrons in our regime, the electron density as a function 
of the current density may be integrated to yield 
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The integration of the nonlinear coupled equations for 
nonrclativistic ions under the self-consistcntly evolving 
potential using Eq. Q] with the boundary condition and 
initial value settings of Mako and Tajima leads to the 
maximum ion (with charge Q) energy as e max = (2a + 
1)QEq, i.e. the maximum carbon energy is determined 
by the maximum electron kinetic energy. In the present 
experimental condition we recognize that the electron ki- 
netic energy is primarily dominated by the laser driven 
collective momentum so that 



(2a + l)Qmc 2 ( 7p - 1) 



(2) 



1 - 1 



which gives rise to (2a + \)Qmc? (^J a\ 

When we insert the exponent of a = 3 as observed 
in our 2D-PIC simulations, we obtain for the case of 
ao = 3.6 a theoretical maximum carbon (proton) energy 
of 59 MeV (9.6 MeV), in the case of a a = 5, 88 MeV 
(14.3 MeV). This is in agreement with the experimental 
value between 36-71 MeV (11-13 MeV). 

Again, for conventional TNSA we do not expect a 
strong dependence of the maximum ion energy on tar- 
get thickness when the target is much smaller than the 
transversal size of the laser focal spot as it is the case in 
the present study. The thickness dependence in TNSA is 
mainly due to transversal spread of the laser accelerated 
electrons on their way through the target 2l|, 127J . 

Instead, the strong increase in e max as it was observed 
in the experiment is related to the growing transmission 
for thinner targets leading to collective electron dynamics 
imposed by the laser field. For thick targets, however, a 
drops towards zero and Eq. 2 reveals the dependence of 
£max on the hot electron energy mc 2 (^Ja 2 L + 1 — 1) [tlli^. 

In summary, ion acceleration from ultra-thin DLC foils 
of thickness 2.9 — 50 nm was studied. A strong depen- 
dence of the resulting maximum energies on target thick- 
ness was observed experimentally, with a pronounced op- 
timum for an initial foil extension of 5.6 nm with a very 
high conversion efficiency from the laser to the kinetic en- 
ergy of the ion bunch. Here, proton energies of 13 MeV 
with a CE of 1.6 % and carbon ions C 6+ of 71 MeV with 



a CE of 10 % were obtained with laser pulse energies as 
low as 0.7 J. Our experimental results are in excellent 
agreement with 2D-PIC simulations. However, we find 
that previously published TNSA scaling laws based on 
analytical models 2^, 21 , 22| fail to interpret our results. 
Therefore, a new semi-analytical model based on [25| was 
developed, showing good agreement with the experimen- 
tal data. 

Demonstrating enhanced ion energies and CE from 
table-top laser systems operating at high repetition rates 
represents a major step towards potential applications 
such as medical or industrial motivated approaches. 
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